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1. Introduction 

Imagine a big birthday cake for Arkady Vainslitein, eacli candle on that 
cake corresponding to one of his outstanding contributions to the modern 
particle theory. I think, a very bright and illuminating candle should then 
mark QCD sum rules. 

The renown papers introducing QCD sum rules [ 1] have been published 
by Shifman, Vainshtein and Zakharov in 1979. The method, known also 
under the nickname of SVZ or ITEP sum rules, very soon became quite 
popular in the particle theory community, especially in Russia. Not only 
experienced theorists, but also many students of that time contributed to 
the development of this field with important results. It was indeed a lot 
of fun to start with an explicit QCD calculation in terms of quark-gluon 
Feynman diagrams and end up estimating dynamical characteristics of real 
hadrons. The flexibility and universality of the sum rule method allowed 
one to go from one interesting problem to another, describing, in the same 
framework, very different hadronic objects, from pions and nucleons to 
charmonium and B mesons. Nowadays, QCD sum rules are still being 
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actively used, providing many important applications and representing an 
important branch on the evolution tree of approximate QCD methods. 

In this short overview, I start, in Sect. 2, from explaining the basic idea 
of sum rules which is rooted in quantum mechanics. After that, in Sect. 3, I 
outline the SVZ sum rule derivation in QCD. Some important applications 
and extensions of the method are listed in Sect. 4. Furthermore, in Sect. 5 I 
demonstrate how QCD sum rules are used to calculate the soft contributions 
to the pion form factor. The light-cone version of sum rules is introduced. 
Many interesting applications of QCD sum rules remain outside this survey, 
some of them can be found in recent reviews [ 2,3]. 

2. SVZ sum rules in quantum mechanics 

To grasp the basic idea of the QCD sum rule method it is sufficient to 
consider a dynamical system much simpler than QCD, that is quantum 
mechanics of a nonrelativistic particle in the potential V{r). The latter has 
to be smooth enough at small distances and confining at large distances. 
The spherically-symmetrical harmonic oscillator V(r) = muP'T^ j1 is a good 
example. Evidently, having defined the potential, one is able to solve the 
problem exactly e.g., by means of the Schrodinger equation, H^pn{x) = 
En'4'ni'r), with the Hamiltonian H = p ^/2m + V{r), obtaining the wave 
functions ijjn{f) and energies En of all eigenstates, n = 0, 1, ... . 

As demonstrated in [ 4] , it is possible to use an alternative procedure 
allowing one to calculate approximately the energy Eq and the wave function 
at zero, V'o(O), of the lowest level. The starting object is the time-evolution 
operator, or the Green's function of the particle G{x2,xi; — ii), taken at 
xi = X2 = and written in terms of the standard spectral representation: 

oo 
n=0 

Performing an analytical continuation of the time variable to imaginary 
values: t2 — t\ ^ —it, one transforms Eq. (1) into a sum over decreasing 
exponents: 

oo 

G(0, 0; -it) = M{t) = |V'n(0)|'e-^"- . (2) 

n=0 

The function M{t) has a dual nature depending on the region of the 
variable r. At small r, the perturbative expansion for M(t) is valid, and 
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it is sufficient to retain a few first terms: 
Mf^'^*(T) = Mf'-'^^T) ^1 - 4:mJrdrV{r)e-'^""''/'' + ©(V^) + ...^ , (3) 

where Mf'^^^{T) = (^)^/^ is the Green's function of the free particle 
motion. Using QCD terminology, we may call the behavior of M(t) at 

small T "asymptotically free" having in mind that it is approximated by a 
universal, interaction-free particle motion. Equating (2) and (3) we obtain 

00 

^|V'n(0)re-^''-~Mf-*(r), (4) 

a typical sum rule which is valid at small r, relating the sum over the bound- 
state contributions to the result of the pcrturbativc expansion. Note that 
the latter includes certain "nonperturbative" or "long-distance" effects too, 
namely the subleading terms containing the interaction potential V. 

At large r one has a completely different picture. In the spectral rep- 
resentation (2) the entire sum over excited levels dies away exponentially 
with respect to the lowest level contribution: 

lim M(r) = IVoWpe-^"^. (5) 

r— >oo 

Thus, at large r one encounters a typical "confinement" regime, because 

the lowest level parameters determining M{t) essentially depend on the 
long-distance dynamics (in this case determined by V{r) at large r). 

An important observation made in [ 4] is that at intermediate values of 
r both descriptions (3) and (5) are approximately valid (see Fig. 1). It is 
therefore possible to retain only the lowest-level contribution in the sum 
rule (4) allowing one to estimate both Eq and |^/'o(0)|, without actually 
solving the Schrodinger equation. 

To further improve the quality of this determination, Mp^'^*{t) can be 
rewritten in a form of the integral 

00 

MP''''\t) = j ff"'''\E)e-^^dE, (6) 


resembling the spectral representation, so that the positive definite function 
pP^^*[E) can be called perturbative spectral density. The integral in Eq. (6) 
is then splitted into two parts, introducing some threshold energy Eth > Eq 
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and the sum over all excited states n > 1 in Eq. (4) is approximated by an 
integral over pP^^' (i?) starting from this threshold: 

J2\M0)fe-''-^ - lf'''\E)e-^-dE. (7) 

n=l / 

The latter equation can be called a "duality" relation having in mind duality 
between the asymptotic-freedom regime and the spectral sum. The integral 
(7) is then subtracted from both sides of Eq. (4) leading to the sum rule 
for the lowest level: 

IV'oWl^e-^"^ = j pP^^\E)e-^^dE. (8) 



Note that one could make use of the sum rule relations similar to Eq. (4) 
in an opposite way. Imagine that the interac;tion potential is unknown but 
we have a possibility to measure, for a set of low levels, their wave functions 
at zero separations and energies experimentally. The sum rule (4) could 
then be used to extract or at least constrain the potential V{r). 

Interestingly, quantum mechanics may also serve as a model for more 
complicated patterns of nonperturbative interactions, in which case sum 
rules have to be treated with care. An example presented in [5] is a po- 




Figurc 1. The analytically continued Green's function M{t) for a particle in the oscilla- 
tor V{r) = moj^r^ /2, normalized to M-I''^'"^{t) and plotted as a function of tui. The exact 
solution (solid) is compared with the perturbative calculation including the first order 
in V correction (dashed) and with the contribution of the lowest bound state (dotted). 
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tential containing two terms: V{r) = Vb[(r/r4)^ + (r/rn)^^] with m <C r4, 
that is, a sharp short-range confining potential combined with a broader 
one. At T — > 0, in the perturbative expansion of M(t) the correction due 
to the second term in the potential is much smaller than the correction 
associated with the first term. However, if one ignores the "nonperturba- 
tive effect" related to the short-distance scale rn, the resulting sum rule 
simply reproduces the lowest level in the potential V{r) = Vb(r/r4)'*. In 
reality, the physical picture is quite difi^erent because it is the ^ r^^ part 
of the potential which mainly determines the formation of bound states. 
The important lesson drawn from this example is: if for some reason a 
short-distance nonpcrturbativc effect is missing and/or ignored, the sum 
rule does not work (or, in other words, duality is violated). 

Interestingly, the sum rule approach in quantum mechanics can be gen- 
eralized to calculate more complicated characteristics such as the ampli- 
tudes of electric-dipole transitions between the lowest S and P levels in 
a given nonrelativistic potential [ 6] . One has to construct a three-point 
correlation function: 



where v = i{Hx — xH) is the quantum-mechanical velocity operator, and 
e- V is the operator corresponding to the dipole radiation of a photon with 
polarization e. The correlator (9) corresponds to the propagation of a 
particle in P wave (below threshold or in imaginary time) from point 2 to 
point 3 where a dipole photon is radiated and tivm further propagation to 
point 1 in 5 wave. Calculating Eq. (9) perturbativcly and matching it to 
the double spectral sum over P and S levels, one gets a sum rule which, 
at intermediate values of the two variables ti,2 is well approximated by the 
contributions of the three lowest El transition amplitudes (IP IS, 2S 



The sum rule approach considered here is, of course not very important 
for quantum mechanics itself, but as we shall see in a moment, serves as a 
very convenient prototype for an analogous method in QCD, in the theory 
where no exact solution is so far available. 

3. SVZ Sum rules in QCD 

We now move from the safe haven of nonrelativistic quantum mechanics 
to QCD, a complicated theory with a rich pattern of quark- gluon and 
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gluon-gluon interactions. At short distances, due to asymptotic freedom 
the theory can stiU be resolved. One considers a quasi-free quark prop- 
agation with calculable pcatiirbative corrections. However, at large dis- 
tances, r ~ 1/Aqc£), the QCD perturbation theory becomes inapplicable 
and the confinement phenomenon takes over, driven by the quark-gluon 
fluctuations in the QCD vacuum. As a result, quarks build coherent bound 
states, hadrons. In general, it is not possible to describe QCD interactions 
with a potential. Nevertheless, qualitatively, the pattern of quark-antiquark 
forces in QCD, with asymptotic freedom at small distances and formation of 
bound states at large distances, is very similar to the quantum-mechanical 
motion in the confining, oscillator-type potential considered in the previous 
section. It is therefore not surprising that sum rules [ 1] analogous to the 
quantum- mechanical ones exist also in QCD. 

The starting object in QCD analogous to the Green's fimction G{0, 0, t) 
is the correlation function describing an evolution of a colorless quark- 
antiquark pair emitted and absorbed by external currents. A "classical 
example" considered in [ 1] is the correlation of two = {u'-f^u — dri^d)/2 
quark currents with the p° meson quantum numbers (isospin 1, = 1~): 

n^.(?) = i J d^xe"i''{0\T{jP{x),jm}\0) ■ (10) 

The dispersion relation (Kallen-Lehmann representation) for this correla- 
tion function contains a sum over all intermediate hadronic states, a direct 
analog of the spectral representation (1) : 

mr — q'' 

h 

Note that, for brevity, I wrote the above relation in a very schematic way, 
including the excited p resonances and the continuum states with p quan- 
tum numbers in one discrete sum. 

The Borel transformation, B{l/(m^j — g^)}— >exp(— m^^/Af^), converts 
the hadronic representation (11) into a sum over decreasing exponents, 

BU,^ = l](0|j^|/i)(/ib-.lO)e— ^-/^^ , (12) 

h 

i.e., the inverse Borel variable 1/M^ plays essentially the same role as the 
auxiliary variable r in the quantum-mechanical case. Another very impor- 
tant virtue of the Borel transformation is that it kills subtraction terms in 
the dispersion relation. 

At large spacelike momentum transfers < Q, Q'^ = —<f ^ J^cd 
(corresponding to large M Aqcd after Borel transformation) the quark- 
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antiquark propagation described by the correlation function (10) is highly 
virtual, the characteristic times/distances being a;o ~ |a;| ~ 1/\/Q^. One 
can then benefit from asymptotic freedom and calculate the correlation 
function in this region porturbatively. The corresponding diagrams up to 
0{as) are depicted in Fig. 2. 

As first realized in [1], there are additional important effects due to 
the interactions with the vacuum quark and gluon fields. The latter have 
typically long-distance (~ ^qcd) scales and, in first approximation, can be 
replaced by static fields, the vacuum condensates. An adequate framework 
to include these effects in the correlation function was developed in a form 
of the Wilson operator prodiict expansion (OPE). The Borel transformed 
answer for the correlation function (10) reads: 

Bn^J'^ = Bnp*+ BCI^MOA^), (i3) 

d=3,4,... 

where the first term on the r.h.s. corresponds to the perturbative diagrams 
in Fig. 2, whereas the sum contains the contributions of vacuum conden- 
sates, ordered by their dimension d. Diagrammatically, these contributions 
are depicted in Fig. 3. The terms with d < 6 contain the vacuum aver- 




(b) (c) (d) 

Figure 2. Diagrams determining the perturbative part of the correlation function (10): 
the free-quark loop (a) and the O(os) corrections (b,c,d). Solid lines denote quarks, 
dashed lines gluons, wavy lines external currents. 

ages of the operators O3 = qq, O4, = G^^G"^'', O5 = gcT^i.(A»/2)G'^^'^g, 
Oe = {qTrq){qT,q), and Of = fabcG''^^G''^''G''^i', where q = u,d,save the 
light-quark fields, G^^ is the gluon field strength tensor, and r,.,s denote 
various combinations of Lorentz and color matrices. Importantly, to com- 
pensate the growing dimension of the operator O^, the Wilson coefficients 
G^^ contain increasing powers of Correspondingly, BC^^^ contain 

powers of 1/M^, making it possible at large to retain in the r.h.s. of 
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Eq. (13) only a few first condensates. Tlius, at ~ 1 GeV^ it is practi- 
cally possible to neglect all operators with d> 6. 

Equating at large the hadronic representation to the result of the 
OPE calculation we obtain the desired sum rule: 

J2iVf')iW)e-"'''^^''^ = BU^T + E BC^,A0\O,\0). (14) 

h d=3,4,-- 

The explicit form of this relation is [ 1] : 

fpe~"^''^'^ + {excited, contimmm p states} 



1 I as{M) \ ^ {mu + md){qq) 



47r2 V 71" y 

1 {^Gl^G-n 1127ra,(g9)2 



^12 M4 81 MR J ' ^^^^ 

where the decay constant of the p meson is defined in the standard way, (p*^ | 
jP I 0) = p I ^fi)ra pt'v^* . In obtaining this relation the four-quark vacuum 
densities are factorized into a product of quark condensates. The quark- 
gluon and three-gluon condensates have very small Wilson coefficients and 
are neglected. The strong coupling a., is taken at the scale M which is the 
characteristic virtuality of the loop diagrams after the Borel transformation. 
A more detailed derivation of this sum rule can be found, e.g. in the review 
[ 3] . The QCD vacuum condensates were recently discussed in [ 7] . 

In full analogy with quantum mechanics, there exists a SVZ region of 
intermediate where the p meson contribution alone saturates the l.h.s. 
of the sum rule (15). To illustrate this statement numerically, in Fig. 4 
the experimentally measured fp (obtained from the p° e+e" width) 
is compared with the same hadronic parameter calculated from Eq. (15) 




(d) (e) (f) 

Figure 3. Diagrams corresponding to the gluon (a,b,c,), quark (d), quark-gluon (e) and 
four-quark (f) condensate contributions to the correlation function (10). The crosses 
denote the vacuum fields. 
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where all contributions of excited and continuum states are neglected. One 
indeed observes a good agreement in the region 1 GeV^. 




Figure 4. The p meson decay constant calculated from the sum rule (15) neglecting all 
excited and continuum states (solid), as a function of the Borcl parameter, in comparison 
with the experimental value (boxes). The dashed curve corresponds to an improved 
calculation , where the sum over excited and continuum states is estimated using quark- 
hadron duality with a threshold Sg = 1.7 GeV^. 



An important step to improve the sum rule (14) is to use the quark- 
hadron duality approximation. The perturbative contribution to the corre- 
lation function (the sum of Fig. 2 diagrams) is represented in the form of a 
dispersion integral and splitted into two parts: 

BIL^r = J f^;j'e-^/^' + j if;:'e-''^' . (16) 

so 

The sum over excited state and continuum contributions in Eq. (14) is ap- 
proximated by the second integral over the perturbative spectral density 
^pert rpj^jg integral is then subtracted from both parts of Eq. (14). Cor- 
respondingly Eq. (15) is modified: the l.h.s. contains only the p term , 
and, on the r.h.s., the perturbative contribution has to be multiplied by 
a factor (1 — e~*"/^^ ). The numerical result obtained from the duality 
improved SVZ sum rule is also shown in Fig. 4. Quark-hadron duality can 
independently be checked for the channels with sufficient experimental in- 
formation on excited hadronic states, such as the J/ij) and p channels. For 
V' resonances one of the first analyses of that type has been done in [ 8] . 

Importantly, not all correlation functions lead to valid QCD sum rules. 
The response of the QCD vacuum to the quarks and gluons "injected" by 
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external currents crucially depends on the quantum numbers and flavour 
content of the current. After all, that is the main reason why hadrons 
are not alike [5]. In certain channels, e.g. for the correlation functions 
of spin zero light-quark currents, specific short-distance nonperturbative 
effects related to instantons are present (the so called "direct instantons" ) [ 
5,9]. These effects remain important even at comparatively large and 
are not accountable in a form of OPE. The subtle duality balance between 
"quasiperturbative" OPE and resonances is destroyed in such cases. Due to 
instantons, it is not possible, for example, to calculate the pion parameters 
using correlators with pscudoscalar uj5d currents. Models based on the 
instanton calculus have to be invoked (sec ;C.g.[ 10]). 

The QCD procedure outlined in this section has indeed many similarities 
with the sum rule derivation in quantum mechanics. To make the analogy 
more transparent, in the following table I put together the main points of 
the two sum rule approaches: in quantum mechanics and in QCD. 



Quantum mechanics 


QCD 


Particle in a smooth confining potential 


quark-antiquark pair in QCD vacuum 


Green's function G{0,0,t) 


Correlation fiuictiou ll/j,y(<:/'^) 


Spectral representation 


Dispersion relation in q'' 


Analytical continuation to — ir 


Borel transformation q'^ — > M'^ 


Perturb, expansion in powers of V 


OPE (Condensate expansion) 


matching the lowest level to M*"^'"* (t) 


matching the lowest hadron to BTl{M^) 


duality of the quasifrce-motion 
and the spectral sum 


quark-hadron duality 


extracting V{r) 
from exp. known spectral sum 


extracting condensates, mu,d,s,c,b 
from exp. known spectral density 


sum rule does not work 
if the short-distance part 
of the potential is ignored 


sum rules do not work 

in the channels 
with direct instantons 


3-point sum rules 
for El-transition amplitudes 


3-point sum rules 
for hadronic matrix elements 



4. Applying and extending the method 
4.1. Baryons 

Following very successful applications of QCD sum rules in the mesonic 
channels [ 1], the next essential step was to extend the method to the 
baryonic sector [11,12]. Correspondingly, the correlators of specially con- 
structed quark currents with baryon quantum numbers were considered. A 
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well known example is the loffe current with the proton quantum numbers: 
J^(x) = eabc{u''^{x)C-i^u\x))^5rd'{x) , (17) 

where a,b,c are color indices and C is the charge conjugation matrix. From 
the QCD sum rule for the correlator (0| JAr(a;) jj^(0)|0) an approximate for- 
mula can be obtained, 

TUN ^ [-(2.0)(27r)2(0|5g|0)(/x = IGeV)]^/^ , (18) 

relating the nucleon mass and the quark condensate density. Thus, QCD 
sum rules unambiguously confirm the fundamental fact that ~ 99% of the 
baryonic mass in the Universe is due to the vacuum condensates. 



4.2. Quark mass determination 

The sum rule relations similar to Eq. (14) are widely used to extract the 
fundamental QCD parameters, not only the condensates themselves but 
also the quark masses. One needs sufficient experimental data on hadronic 
parameters in a given channel (masses and decay constants of ground and 
excited states, experimentally fitted ansatze for continuum states) in order 
to saturate the hadronic part of the sum rule. 

The ratios of the light {u, d, s) quark masses are predicted from the QCD 
chiral perturbation theory [ 13]: 

— = 0.553 ± 0.043, — = 18.9 ±0.8, 2— = 24.4 ± 1.5 , (19) 

rud ma m„ + 

(there is also a more recent estimate mu/md = 0.46 ± 0.09 [ 14]). QCD 
sum rules offer a unique opportunity to estimate the individual masses of 
u, d, s quarks. To illustrate the continuous efforts in this direction, let me 
mention one recent determination of the strange quark mass [15], based 
on the correlation function of the derivatives of the strangeness-changing 
vector current = S'y^q, q = u,d: 

W{q) = i J d^xe"'^0\T{d^j^^{x)d^j^^O)}\0) . (20) 

The OPE answer for 11* is proportional to (wg — m^)^ ~ turning this 
correlator into a very convenient object for the mg extraction. Furthermore, 
the recent progress in the multiloop QCD calculations allows to roach the 
0{a^) accuracy in the perturbative part of 11''. An updated analysis of 
kaon S wave scattering on tt, rj, rj' is used to reproduce the hadronic spectral 
density. The sum rule yields [ 15] for the running mass in the MS scheme: 
ms(2GeV) = 99 ± 16 MeV, in a good agreement with the recent lattice 
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QCD estimates. Using the ratios (19) one obtains m„(2GcV) — 2.9 ± 0.6 
MeV and md(2GeV) = 5.2 ± 0.9 MeV. The earlier work on predicting the 
light-quark masses from QCD sum rules is summarized in [ 3,16]. 

The charmed quark mass determination was one of the first successful 
applications of the QCD sum rule approach [ 17,1]. The correlation function 
of two CTfpiC currents (in other words, the charm contribution to the pho- 
ton polarization operator) was matched to its hadronic dispersion relation, 
where the imaginary part is simply proportional to the e~ — » charm 
cross section including ^ resonances and the open charm continuum. The 
lowest power moments of this sum rule at = are well suited for rric de- 
termination because nonpertrbative effects are extremely small. Replacing 
c ^ b, ip ^ T and open charm by open beauty one obtains analogous sum 
rule relations for the b quark [ 18]. In recent years the mainstream devel- 
opment in the heavy quark mass determination went in another direction, 
employing the higher moments which are less sensitive to the experimental 
input above the open flavour threshold. These moments, however, demand 
careful treatment of Coulomb corrections [ 19] which is only possible in the 
nonrclativistic QCD (the current status of this field is reviewed in [20]). Re- 
cent precise measurements of the e+e~ hadrons cross section on one side 
and a substantial progress in the calculation of perturbative diagrams on 
the other side, allowed to reanalyze with a higher precision the low moments 
of the original SVZ sum rules for quarkonia with the following results [ 21] 
for the MS masses: TOc(mc) = 1.304 ± 0.027 GeV, TOf,(mb) = 4.209 ± 0.05 
GeV. Another subset of charmonium sum rules (higher moments at fixed 
large q"^ < 0) was employed in [ 22], with a prediction for rric in agreement 
with the above. 

4.3. Calculation of the B meson decay constant 

Having determined the condensates and quark masses from a set of exper- 
imentally proven QCD sum rules for light-quark and quarkonium systems 
one has an exciting possibility to predict the unknown hadronic characteris- 
tics of B meson. In the amplitudes of exclusive weak B decays the hadronic 
matrix elements are multiplied by poorly known CKM parameters, such as 
Vub- QCD sum rule calculations may therefore provide a useful hadronic 
input for extraction of CKM parameters from data on exclusive B decays. 
Importantly, the theoretical accuracy of the sum rule determination can be 
estimated by varying the input within allowed intervals. 

One of the most important parameters involved in B physics is the B 
meson decay constant fs defined via the matrix element {0\uij5b\B) . The 
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calculation of fs using QCD sum rules has a long history, a detailed review 
and relevant references can be found, e,g. in [3,23], I only mention the 
very first papers [ 17,24]. One usually employs the SVZ sum rule for the 
two-point correlator of 6175(7 currents. I will not write down this sum rule 
explicitly. It looks very similar to the one for fp discussed in sect. 3, in a 
sense that the sum rule contains a (duality subtracted) perturbative part 
and condensate terms. The expressions for the Wilson coefficients are in 
this case much more complicated, especially the radiative corrections to the 
heavy-light loop diagrams. The recent essential update of the sum rule for 
fs is worked out in [ 25] taking into account the 0{a^) corrections to the 
heavy-light loop recently calculated in [ 26] and treating the b quark mass 
in MS scheme. The result is (for mb{mb) = 4.21 ±0.05 GeV): /b = 210 ±19 
MeV and fs^ = 244 ± 21 MeV, in a good agreement with the most recent 
lattice QCD determination (including dynamical sea-quark effects). I think, 
the example of Jb determination demonstrates that QCD sum rules indeed 
provide a reliable analytical tool for the hadronic B physics. 



4.4. Hadronic amplitudes 

To complete this short survey of QCD sum rule applications, it is important 
to mention that this method allows to calculate various hadronic amplitudes 
involving more than one hadron. Let me consider, as a generic example, a 
calculation of the hadronic matrix element {hf{p + q)\j\hi{p)) of a certain 
quark current j with a momentum transfer q. The convenient starting 
object in this case is the three-point correlation function depending on two 
independent 4-momenta: 

Tf,{p, q) = {if I d^x d*y e-'(P-+''yH0\T{jf{0)j{y)Mx)}\0) . (21) 

As a next step, one writes down a double dispersion relation, in the variables 
and {p + qf at fixed q^, expressed in a form similar to Eq. (11): 



TfM 5) = E E /°l'^"';^^^'^lt^^2''^'''t + subtractions , (22) 

where the double sum includes all possible transitions between the states 
with hi and hf quantum numbers. Two independent Borel transforma- 
tions in and {p + qf applied to Eq. (22) enhance the ground-state term 
containing the desired matrix element and allow to get rid of subtraction 
terms: 

B^B^Tfi = J2 E(0|i/|/*/)('^/lil'^')('^'li'|0)e"'"'^^^''^'"'"^*^''' , (23) 
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P P+q 



(a) (b) (c) 

Figure 5. Contributions to the 3-point correlation function (21): (a) perturbative, ze- 
roth order in as; (b)-(c) some nonpcrturbative corrections. 

where Mi and M2 are the Borel variables corresponding to and (p + g)^ , 
respectively. On the other hand, the correlator (21) can be computed, in 
terms of perturbative and condensate contributions: 

•' d=3,4,.. 

In the above, the perturbative contribution calculated from the diagram 
in Fig. 5a is represented in a convenient form of double spectral repre- 
sentation with a spectral density p^jY* ■ The Wilson coefficients are 
calculated from the diagrams exemplified in Fig. 5b, c. Importantly, one 
does not introduce new parameters/inputs in this calculation, benefiting 
from the universality of quark/gluon condensates. The above expansion 
is valid at large spacelike external momenta: \ {p + q)'^\ 3> ^qcD' f^'^ 
from the hadronic thresholds in the corresponding channels. Accordingly, 
the squared momentum transfer is also kept large, = —q^ 3> ^qcd^ 
at least for the currents containing light quarks Equating the hadronic 
dispersion relation to the OPE result at large , M| and invoking quark- 
hadron duality one obtains the sum rule for the matrix element: 

+ Y,BiB2C%{0\0,\0), (25) 

d 

where fi = {hi\ji\0), ff = {0\jf\hf) are the decay constants of the initial 
and final hadrons. The latter are calculable from two-point sum rules, or 

^Hadronic matrix elements at = 0, e.g., the nucleon magnetic moment can also 
be calculated within QCD sum rule approach, using the external (background) field 
technique [ 27], with additional vacuum condensates induced by external, non-QCD 
fields. 
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simply known from experiment. In the above, R{so, Sq) is the quark-hadron 
duality domain in the s,s' plane, So,Sq are the corresponding thresholds. 
Using 3-point correlators, the sum rules for charmonium radiative transi- 
tions have been derived in [ 28] . Another important application [ 29] is the 
pion e.m. form factor discussed in more detail in the next section. 



5. QCD sum rules and the pion form factor 

One of the celebrated study objects in hadronic physics is the pion 
electromagnetic form factor F^(g^) determining the pion matrix element 
(7r(p + q)\jtrHp)) = F^i<f){'^P + Q)ti of the quark e.m. current fj^ = 

At very large values of the spacclikc momentum transfer = —q^ oo 
the form factor is determined by the perturbative QCD factorization [ 30]: 

2 



FAQ') '^-^'^^ 



9g2 



1 







(26) 



obtained by the convolution of distribution amplitudes (DA) ^Triu) of the 
initial and final pions (see the definition below) with the 0{as) quark-gluon 
hard-scattering amplitude. At finite Q^, the major problem is to estimate 
the "soft", 0{a^g/Q'^) part of this form factor. It corresponds to an overlap 
of end-point configurations of the quark-parton momenta in the initial and 
final pions, so that the large momentum is transferred without a hard gluon 
exchange (the so called Feynman mechanism). 

The first model-independent estimate of the soft contribution to the 
pion form factor was provided by QCD sum rules [ 29] . The three-point 
correlator (21) was used, with j, j, and jf replaced by j^™, j^s and jj^g, 
respectively, where j^n = u'yulbd is the axial-vector current generating the 
pion state from the vacuum: |0) = —if-nPu- The calculation based 

on OPE and condensates is valid at sufficiently large Q^, practically at ~ 
1 GeV^. The resulting sum rule for the form factor written in the form (25) 
has a rather compact expression: 



s + s 



a- „^ . 2087r _ , 2Q2 



where the perturbative spectral density is 

pf-*(s,/,02) = ^J_[3A(a + g2)(a + 2Q2)_A2-5Q2(^ + Q2^3]^(28) 
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with A = (c + (5^)^ — 4ss' and a ~ s + s' . In Eq. (27) the condensates up 
to d = 6 are included, m„ = 0, Mi — M2 = M and sq = Sq = ~ 0.7 
GeV^. The duaUty threshold is inferred from the two-point sum rule for 
the axial- vector channel [1]. At = 1 3 GcV^, the form factor pre- 
dicted from the sum rule agrees with the experimental data. E.g., compare 
F-TTiQ"^ = IGeV^) ~ 0.3 predicted from Eq. (27) with the most accurate 
CEBAF data [ 31] shown in Fig. 7 below. The good agreement indicates 
that the soft mechanism is the most important one in this region and that 
the 0{aa) hard scattering effect which should dominate at infinitely large 
Q"^ is still a small correction. (The latter corresponds to the gluon ex- 
changes added to the diagram of Fig. 5a). At largo the pcrturbativc 
part of the sum rule (27) has a ~ behavior, in full accordance with our 
expectation for the soft, end-point contribution to the form factor. How- 
ever, the condensate contributions to i^7r(Q^) arc cither Q^-independent or 
grow ^ Q^/M^. A careful look at one of the relevant diagrams in Fig. 
5c reveals the reason of this anomalous behavior. Using local (static field) 
condensate approximation, one implicitly neglects the momenta of vacuum 
quark/gluon fields. The external large momentum p is carried by a single 
quark, which, after the photon absorption, propagates with the momentum 
p + q, so that the contribution of this diagram is independent. Therefore, 
the truncated local condensate expansion is not an adequate approximation 
to reproduce the large behavior of the pion form factor. 

A possibility to calculate Ft^{Q'^) including both soft and hard scattering 
effects at large [ 32,33], is provided by the light-cone sum rule (LCSR) 
approach [ 34] combining the elements of the theory of hard exclusive pro- 
cesses [ 30] with the SVZ procedure. 

One starts with introducing a vacuum-to-pion correlation function 



where one of the pions is put on-shell, = m^, and the second one is 
replaced by the generating current j^^. For this correlator a dispersion 
relation is written, in full analogy with Eq. (11): 




(29) 




(30) 



The lowest pion-state term {h = n) in the hadronic sum, 




(31) 
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p + g 9 



(a) (b) (c) 

Figure 6. Light-cone expansion of the correlation function (29). 

At large spacelike momenta, + ^ ^qcdj the correlation 

function (29) is dominated by small values of the space-time interval x^, 
allowing one to expand the product of two currents around the light-cone 

= 0. The leading-order contribution is obtained from the diagram in 
Fig. 6a and consists of two parts: (1) the short-distance amplitude involving 
the virtual quark propagating between the points x and 0, and (2) the 
vacuum-to-pion matrix element of a nonlocal quark- antiquark operator, 
{Q\u{x)^ p^c^d{Q)\-K) . This matrix element contains long-distance effects and 
is therefore not directly calculable. On the other hand, being expanded 
near = it can be resolved in terms of universal distribution functions: 

(0|M(a:)7^75d(0)|7r(p)) = ~ip^,f^ [ due^^^-"^ {(p^{u, n) + x^gi{u, fi)) 

Jo 

+ fJx^-^) [ due-'^P-^g2iu,ii) + ...,{32) 
V P ' X / Jo 

where the terms up to 0{x^) are shown explicitly. This expansion contains 
normalized light-cone distribution amplitudes (DA): (pT^{u,fi), gi{u,fi), 
(72 (u, A*), ■■• 7 and the scale ^ reflects the logarithmic dependence on x^. 
Importantly, the power moments of DA, e.g., Mn{fJi) = Jq du w'V7r(w, /x), 
are related to the vacuum-to-pion matrix elements of local quark-antiquark 
operators with a definite twist (dimension minus Lorentz spin). For that 
reason, is called twist 2 DA, and, correspondingly, 171 and g2 are of twist 
4. Thus, in the light-cone OPE one deals with a completely different pat- 
tern of long-distance effects, as compared with the local OPE considered 
in Sect. 3. Instead of a set of universal vacuum condensates, there is a set 
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of DA for a given light meson, each of DA representing a series of matrix 
elements. 

Actually, the twist 2 DA t^^ was originally introduced in the QCD anal- 
ysis of hard exclusive hadronic processes [ 30], see e.g., Eq. (26). Some of 
its properties are well understood, in particular, the following expansion 
can be written: 

M^, = 6w(l -u)(^l + a2n{tJ-)Cti\2u - 1)^ , (33) 

based on the approximate conformal symmetry of QCD with light quarks. 
In the above, C2n are Gegenbauer polynomials and the coefficients a2„(/x) 
determine the deviation of (/?^ from its asymptotic form 6u{l — u). Due to 
the perturbative evolution, a2,i(/i) arc logarithmically suppressed at large 
IJ,. The low-scale values of a2n (and of similar coefficients for other DA) 
have to be considered a nonperturbative input. 

The correlation function (29) calculated from the light-cone OPE rep- 
resents a convolution of the pion DA and short-distance (hard scattering) 
amplitudes: 

1 

F,.{p, q) = 2^Up,p. j du^-^^^^^—^^ +.... (34) 



For simplicity, only the leading order, twist 2 contribution with the rele- 
vant kincmatical structure is shown. The ellipses denote the 0{as) cor- 
rections (one of the diagrams is presented in Fig. 6c) and the higher twist 
contributions suppressed by powers of the denominator. Physically, the 
higher-twist corrections take into account the transverse momentum of the 
quark- antiquark state (e.g., the twist 4 terms in the expansion (32)) and 
the contributions of higher Fock states in the pion wave function (such as 
the quark-antiquark-gluon DA contributing via the diagram in Fig. 6b.). 
These two effects are related via QCD equations of motion. More details 
on the pion DA of higher twists can be foimd in [ 35] . The factorization 
scale /X in Eq. (34) effectively separates the large virtualities (> /x^) in the 
hard scattering amplitude from the small ones (< /it^) in the pion DA. 

Equating the dispersion relation (30) to the OPE result (34) at large 
spacelikc (p+q)^: one extracts the form factor Ft^{Q'^) applying the standard 
elements of the QCD sum rule technique, the Borel transformation in the 
variable {p + q)^ and the quark- hadron duality. The latter reduces to a 
simple replacement of the lower limit in the u-integration in Eq. (34), 
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Figure 7. The pion cm. form factor calculated from LCSR in comparison with the 
CEBAF data shown with points. The solid line corresponds to the asymptotic pion DA, 
dashed lines indicate the estimated overall theoretical uncertainty ; the dash-dotted line 
is calculated with the CZ model of the pion DA. 

Q2/(sJ + rpjjg resulting sum rule [ 32,33] is: 
1 

QV(sJ+Q^) 

where the leading-order twist 2 part is shown explicitly. The behavior 
of Eq. (35) corresponds to the soft end-point mechanism, provided that in 
the oo limit the integration region shrinks to the point u = 1. The 

0{as) part of this sum rule was calculated in [33] and is indeed small 
at 0^ _ 1 GeV^. Importantly, in Fi"^^(Q2) 

one recovers the ~ 
asymptotic term corresponding to the hard perturbative mechanism, with 
a coefBcient which, in the adopted approximation, coincides with the one in 
Eq. (26). The higher twist contributions to the sum rule (35) manifest the 
same ^ l/Q^ behavior as the leading twist. Altogether, we seem to achieve 
the goal. The pion form factor obtained from LCSR contains both the hard- 
scattering and soft (end-point) contributions, with a proper asymptotic 
behavior at large Q^. 

The updated LCSR prediction for Ft^{Q'^) [36] is shown in Fig. 7. 
One important practical use of this result is to estimate/constrain the 
nonasymptotic coefficients a-2n by fitting the sum rule (35) to the ex- 
perimental data on the pion form factor. However, currently there are 
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no sufficient data at >1 GeV^ to constrain complicated patterns of 
nonasymptotic coefficients. Considering simple ones, one finds that, e.g. 
the asymptotic DA (pn{u) = Qu{l — u) is not excluded, whereas the 
CZ -model [ 37] seems to be disfavored by data. Assuming that only 
a2 ^ and neglecting all other coefficients yields [ 36] the following range 
02(1 GeV) = 0.24 ± 0.14 ± 0.08, where the first error reflects the estimated 
theoretical uncertainty and the second one corresponds to the experimental 
errors. 

Other LCSR applications to the physics of hard exclusive processes in- 
clude the 7*77r form factor [38], the kaon e.m. form factor [36], and the 
first attempt to calculate the nucleon form factors [ 39] . 

Furthermore, an important task of LCSR is to provide B physics with 
various heavy-to-light hadronic matrix elements. In particular, the sum rule 
for the B TT form factor can be obtained from a correlation function very 
similar to the one depicted in Fig. 6, if one replaces the virtual light quark by 
a b quark. The calculable short-distance part will then change considerably, 
but the long-distance part remains essentially the same, determined by 
the set of pion DA. The sum rule predictions for the B ^ tt [ 40] and 
B ^ p [41] form factors are already used to extract \Vub\ from the widths 
of B ^ 'k{p)Iui decays. One can also employ LCSR to estimate the hadronic 
amplitudes for B tttt and similar decays beyond factorization [ 42] . A 
summary of the sum rule applications to the heavy flavour physics can be 
found in [ 43] . 

6. Conclusion 

For more than twenty years, QCD(SVZ) sum rules serve as a virtual lab- 
oratory for studying the transition from short- to long-distance QCD. The 
practical use of this analytical method is twofold. On one hand, using sum 
rules for experimentally known hadronic quantities, QCD parameters such 
as quark masses are extracted. On the other hand, the sum rules are em- 
ployed to predict unknown hadronic parameters, for example /s, with a 
controllable accuracy. Finally, the example of the pion form factor calcula- 
tion demonstrates that the light-cone version of QCD sum rules has a large 
potential in describing exclusive hadronic transitions. 
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